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Abstract 



We consider the properties of the ground state of bottomium. The T 
mass is evaluated to two loops, and including leading higher order [O(aslogas)] 
and m^/ml corrections. This allows us to present updated values for the pole 
mass and MS mass of the b quark: mi, — 5022 ± 58 MeV, for the pole mass, and 
fhb{fhb) = 4286 ±36 MeV for the MS one. The value for the MS mass is accurate 
including and 0{a'l) corrections and leading orders in the ratio mi/m^. We then 
consider the wave function for the ground state of bh, which is calculated to two 
loops in the nonrelativistic approximation. Taking into account the evaluation of 
the matching coefficients by Beneke and Signer one can calculate, in principle, the 
width for the decay T — > e^e~ to order . Unfortunately, given the size of the 
corrections it is impossible to produce reliable numbers. The situation is slightly 
better for the ground state of toponium, where a decay width into e^e~ of 11 - 
14 keV is predicted. 
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-THE GROUND STATE OF BOTTOMIUM TO TWO LOOPS AND HIGHER- 

1. Introduction 

In this paper we profit from tlie existence of recently obtained results on heavy quarkonium to give an 

accurate picture of the ground state of bottomium. In particular, we have the two loop relation between 
pole and MS massI^J; the O(aglogas) and 0{ml/ml) corrections to quarkonium potential I^' ; and the two 
loop expression for the decay T — > e+e~ in terms of the nonrelativistic ground state wave functionl^l, which 
last we evaluate here also to two loops. 

All this allows us to give a determination of the pole mass of the b quark, 

nib = 5022 ± 58 MeV 

correct up to and including and 0{al logag) corrections as well as leading nonperturbative and 0{ml/ml) 
ones; to find a value for the MS mass, 

rfibiffib) = 4286 ± 36 MeV 

correct to order and leading nonperturbative and 0(m^/m^); and to present an evaluation of the width 

r{r e+e") 

accurate to order and including leading nonperturbative contributions. The determination of the pole 
mass turns out to be very stable, in particular when comparing it with older determinations!'*"''! The MS 
mass is less stable, but on the other hand is in good agreement with recent NNLO (next-to-next to leading 
order) evaluations based on QCD sum rules'^l. Finally, the values found for the width, while in rough 
agreement with experiment, show an excessive dependence both on the renormalization point chosen to 
calculate and on the degree of accuracy of the calculation: both due to the very large size of one and two 
loop corrections. 

The plan of this paper is as follows. In Sect. 2 we present the hamiltonian used to obtain the 
results. This is used in Sect. 3 to evaluate pole mass and MS mass for the b quark, and to comment on the 
convergence properties of the series. The decay rates to e+e~ of the T are discussed in Sect. 4. The article 
is finished with an Appendix where some of the evaluations, which may be of interest for other problems, 
are presented in detail. 



2. The nonrelativistic hamiltonian 

In the static limit the bb interaction may be described by a potential. The ensuing hamiltonian may be 
written as _ 

H = H^°^+Hi (2.1) 

where _ 

m^)=2m+^A-^^^; (2.2a) 
m r 

we have put together all coulomb-like pieces of the interaction so that 

,2\ 



5,(/z2)=a,(M^)|l + cW^^+c(2) 



am 



.dll.2] 



ai = 



7E(a,/3o + ^) + g+7i)f +a2 
31Ca - 20TFnf 



36 



- + f Cs] C^T^nf - [f - I6C3] C^T^Uf + ^T|nJ.} 



13.2. 



(2.2b) 



(2.2c) 



- 1 - 



-F. J. YNDURAIN- 



The remaining piece of the interaction Hi is, in the non-relativistic (static) limit, 

^NR ^ yiL.l) ^ yiL,2) ^ y{LL)^ ^2.3a) 

TT 2 r ' 



Cf/Sq"? log ''M 



475-2 J. 

Note that the mass that appears in above formulas is the pole mass, but all other renormalization is carried 
in the MS scheme. For more details on this see rcfs. 4, 5. 

Eqs.(2.2) give what wc will need to calculate the wave function, when taking into account also the 
results of Beneke et al.I^J; but for the spectrum we can do better. First, we have to add one loop velocity 
corrections. For the spin-independent part of the spectrum this is 

^vel.dep. - 2to^2 ' t^-^a) 

bi = ^{Cp — ^Ca) and a spin-dependent piece 

Fhf = ^^^^^sis + l)5{r) (2.4b) 

has to be added for vector states like the T. Then we have the corrections yielding contributions of order 
a^logois, which are the leading 0{al) corrections in the limit of small a^. They have been evaluated in 
ref. 2 and they lead to a correction on the mass of the v of 

S[al\ogas]Ew = -m[CF + |CA]C|.Q!f (log as)/'7r, (2.5) 

for fj, = 2/a = mCFCts- 

The influence of the nonzero mass of the c quark, the only worth considering, will be evaluated 
now. To leading order it only contributes to the bb potential through a c-quark loop in the gluon exchange 
diagram. The momentum space potential generated by a nonzero mass quark through this mechanism is 
then, in the nonrelativistic limit, 

, -ggiggi log <^ '')-^ . (2^6) 

We expand in powers of m^/k^. The zeroth term is already included in (5). The first order correction is'^l 

OcmassV' = • (2.7) 

In x-space, 

S.^^.V = ^f^^r. (2.8) 

TT 

If we are interested in a calculation correct to order and log ag for the bound state energies, 
or order for the decay rate, all these terms have to be treated as first order perturbations except 
that should be computed to second order. ^ Of course the hamiltonian H^^^ can, and will, be solved exactly. 

3. The b quark mass 



^ In fact, this is the most difficult part of the calculation. 
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3.1. Pole mass 

Here ag is to be calculated to three loops: 



An ( /3ilogL , l3flog'L-pl\ogL + p2po-PU 
"4W = — il + ^472 J (-^-i) 



with 



r — l/^rr ^ . _ 1 1 2„ fj _ ino 38„ o _ 2857 5033 ^ j. 325„2 
yl^! P0 = ll-3«/, Pi -WZ- -^Uf, P2 - ^ T8~^f + ~S4^f- 

We take as input parameters 

A{nf = 4, three loops) = 0.283 ± 0.035 GeV [ a;s(M|) ~ 0.117 ± 0.024] 

(ref. 7) and for the gluon condensate, very poorly known, the value (ttsG^) — 0.06 ± 0.02 GcV* . From the 
mass of the T particle we have a very precise determination for the pole mass of the b quark. This determina- 
tion is correct to order af and including leading 0{m^/ml) and leading nonperturbative corrections as well 
as the al corrections proportional to log ; the details of it will be given below. With the renormalization 
point = nibCFCes we have, 

rUb =5022 ± 43 (A) t 5 i{asG^})ll] (vary fj,"^ by 25%) ± 38 (other th. uncert.) 
=5022 ± 58 McV . 

Here we append (A) to the error induced by that of A, and likewise {{asG^)) tags the error due to that of 
the condensate. The error labeled (other th. uncert.) includes also the error evaluated in ref. 8; the rest is 
as in ref. 5. 

We collect in the table the determinations of the b quark mass based on spectroscopy, to increasing 
accuracy. The stability of the numerical values of the pole mass is remarkable: the pole masses all lie within 
each other error bars. The MS ones (that will be discussed in greater detail in next subsection) show more 
spread. It may be noted that the three loop value for m(,(m^), 4286 ± 36, agrees perfectly with the Beneke- 
Signer three loop result'^l based on sum rules, mb{fhl) = 4260 ± 100 MeV. This agreement, however, is less 
satisfactory than it would appear at first sight because of the large size of the 0{al) corrections. We will 
discuss this more in next subsection. 



Reference 


rrife (pole) 


mb{ml) 


mc(pole) 


mdml) 


TY 

PY 
Here 


4971 ± 72 

5065 ± 60 
5022 ± 58 


4401^35 

4455t^^ 
4286 ± 36 


1585 ±20 (*) 
1866t?^^ 


1321 ±30 (*) 
1542^1^3 



Table 1. b and c quark masses. (*) Systematic errors not included. 



TY: Titard and YndurainW . 0{al) plus 0{al)v, O(v^) for m; 0(a?) for rfi. Rescaled for A{nf = 4) = 283 MeV. 
PY: Pineda and Yndurain'^l Full 0{ai) for m; O(a^) for ffi. Rescaled for A{nf = 4) = 283 MeV. 
Here: This calculation, ©(a^), 0{asrn^ /nil) '^'id 0{allogas) for rn; 0{a'l) and 0{a^m^/ml) for m. Values 
not given for the c quark, as the higher order terms are as large as the leading ones. 



We finally remark that the values of nib quoted e.g. in the Table 1 were not obtained solving Eq. (8), 
but solving exactly the coulombic part of the interaction, and perturbing the result (see refs. 4, 5 for details). 
We also note that, in the determinations of nib, the new pieces, 0{asml/nil) and ©(a^logag), have been 
evaluated to first order; in particular, we have included the corresponding shifts in the central values, not in 
the errors. If we included these, the errors would decrease by some 7%. 

3.2. mft — fhb{fhb) connection 

The connection of the pole mass with the MS mass has been known for some time to one and two loops: 
very recently, a three loop evaluation has been completed. Coupling this with the pole mass evaluations, we 
now have an order result for the MS mass. We review here briefly this. 
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Write, for a heavy quark, 

m{m) = m/{l + Si+S2 + S3-\ }; (3.3a) 

m here denotes the pole mass, and fh is the MS one. One has 

6, = C,^, <5, = c.(^)', ^3 = C3(^)'. (3.3b) 

The coefficient C2 has been evaluated by Gray et al.^^l, and reads 

C2 = -K + 2Cf, (3.3c) 



(3.3d) 



16.11- 1.04 n/; Z\(p) = | [iyr^ - fp" 



3 L8 ' 4' 

rrii are the (pole) masses of the quarks strictly hghter than m, and n/ is the number of these. For the b 
quark case, Uf = A and only the c quark mass has to be considered; we will take mc = 1.8 GeV (see Table 1 
above) for the calculations. 

The coefficient C3 was recently calculated by Mclnikov and van Ritbergenl^l . Neglecting now the m,, 

C3 ~ 190.389 - 26.6551n/ + 0.652694n^. (3.3e) 

For the b, c quarks, with as given before, 





=0.090, 


Siic) 


=0.137, 


S2{b) 


=0.045, 




=0.108, 


Ssib) 


=0.029; 




=0.125. 



(3.4) 

From these values we conclude that, for the c quark, the series has started to diverge at second 
order, and it certainly diverges at order a^. For the b quark the series is at the edge of convergence for the 
al contribution. 

Using the three loop relation (3.3) of the pole mass to the MS mass we then find, from the results 
in Subsect. 3.1 for the pole mass, the value 

fhb{fhb) = 4284 ± 7 (yl) T 5 ((a^G^)) ± 35 (other th. uncert.) = 4284 ± 36 MeV . (3.5) 

This is the value incorporated in Table 1. The slight dependence of to on yl when evaluated in this way was 
already noted in rcf. 4. 

There is another way of obtaining fh, which is to express directly the mass of the T in terms of it, 
using Eq. (3.3) and the order formula for the T mass in terms of the pole mass Eq. (3.3a). One finds, for 
n/ = 4, and neglecting to^/to^. 



,_,.3 (3-6a) 
a., (m) 



M(r) =2m(TO)(l + + 7.559 ( ^ 

66.769+ 18.277 (log Cf+ log as(TO)) ( ^^^^ ) }• 



+ 

(One could add the leading nonperturbative contributions to (3.6a) a la Leutwyler-Voloshin in the standard 
way). This method has been at times advertised as improving the convergence, allegedly because the MS 

mass docs not suffer from nearby renormalon singularities. But a close look to (3.6a) docs not seem to bear 
this out. To an acceptable 0{a'^) error we can replace log(as(TO)) by log(as(M(T/2)) above. With A as 
before (3.6a) then becomes 

M(r) = 2TOb(TOf,) 1 1 + Cf^^ + 7.559 ' + 43.502 (^)' 1 • (3.6b) 
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This does not look particularly convergent, and is certainly not an improvement over the expression using 
the pole mass where one has, for the choice /i = CpmhOis and still neglecting the masses of quarks lighter 
than the b, 

M{T) = 2m, 1 1 - 2.193 ' - 24.725 ' - 458.28 ' + 897.93 [loga.] (^) ' | . 

(3.7) 

To order three, (3.7) is actually better^ than (3.6b). What is more, logarithmic terms appear in (3.6) at 
order a^, while for the pole mass expression they first show up at order a^. Finally, the direct formula for 
M(T) in terms of the MS mass presents the extra difficulty that the nonperturbative contribution becomes 
larger than than what one has for the expression in terms of the pole mass (^ 80 against ^ 9 McV), because 
of the definition of the renormalization point. With the purely perturbative expression (3.6) plus leading 
nonperturbative (gluon condensate) correction one finds the value fhb{fhb) = 4167 MeV, rather low. 



4. The Decay Rates T 



and T 



For bound state calculations it is convenient to solve iJ*^"^ exactly, and treat Hi as a pertiirbation; but, to 
evaluate the wave function it is preferable to work in the following somewhat different manner. Define the 
quantities 

2 2r 



mCFas{ii) '' 



P= —■ 

a 



One may rewrite Hi as 
Hic 



Hi = Hic + h[',^ + H^^l + Hill; 



i asin) , 
2P0 r 



CL + 



7r2 J 



r(2) 

'iz, 



2 CplSoaliij) \ogp 
a 



2tt 



P 



CpalilJ-) logp 



(4.1) 



Hill 



2 CpPiai logp 
a 



4-772 p 

Here cl = aiPo + g/3i + 7e/3o/2- 

One can put H^^^ and Hic together, H^^^ = H^^^ + Hic, and solve the corresponding coulombic 
Schrodinger equation exactly by just replacing, in the ordinary solution of the coulombic problem. 



(4.2) 



{mCp) 



3/2 



V2 



the coupling according to 

ias(/«^) 



1+c' 



(1)! 



+ c 



(2)! 



O^sip) 



Cl 



-rCpnias/'^ 



13^ log fxa/2\ aliiJ.) 



The remaining -ffj^^ +Hill can then be evaluated as first order perturbations; only h[]^ has to be calculated 
to first and second order and, for this second order, the fact that one perturbs on H^°^ has also to be taken 
into account. 



^ The convergence of Eq. (3.7) is still improved if one solves exactly the purely coulombic part of the static potential, 
as was done in Subsect. 3.1. For example, the 0{at) term would become — 232.12(as/7r)''. 
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4.1 Calculation of perturbations 

The NNLO (two loop) hard part of the radiative correction to the leptonic decay of quarkonium, say 
has been evaluated some time ago^l . One can then write 



wherel^l 



r(r ^ e+e") 



QbO^QED 



M{T) 



(1) aJjJ^ (2) Q:Km) 
"hard °v.»-H o 



hard ^2 



-Rio(O) 



(4.3a) 



5(1) 

'hard 

5(2) 
-'hard 



2Cf, 



1 1 

6 — 



79 
36 



log 2 



23 1 /■ 



+CfCa 



151 _ 13 
72 4 



Cs 



+ %CFTpnf + CfTf (-§77^ + f ) 
Cf /3o log m^Z/i^ 



+ 



(4.3b) 



and -Rio(O) is the static (radial) wave function, i.e., evaluated neglecting terms of relative order 1/m in the 
interaction. We will write 

i?,o(0) = |l + + 4^ + ^NP} R'Sm (4.3c) 

i?^o''(0) is as in (4.2) and the (5wf are to be evaluated with H in (2.1) (or (4.1)). Finally, the (leading) 
nonperturbative piece ^np may be obtained in terms of the gluon condensate {agG"^) as in refs. 4, 9 so that 

2 9687r(Q!^G2^ 



3NP 



425m4(CFas)6 " 



Define then, with self-explanatory notation. 



(4.4) 



(4.5) 



Ac is calculated trivially, as it is equivalent to a modification of the coulombic potential in iJ^"^. The 
remaining 6s are easily evaluated with the formulas of the Appendix. We find. 



Ac=l 



Cl 



/3olog/ia/2 



L,2 



Hi 

3/3oas(/x) 



ci + 



/?ologAta/2' 



C2 + 2cl log — + 



a/i , /3glog^a/i/2 



(4.6a) 



47r 



7e + 



+ 3/3o 7e + 



(1) /3ologaM/2 \ q:^(m) 



3 7E + 



TT^ - 6 



2c^+/?o2log^ 



as(M) 
47r2 ' 



(4.6b) 



ii'l=/?o^C3-(l-7E)(l + y 



^21 



7e - 27E + — 
6 



«|0i) 

47r ' 



where A^j^ ^ above actually contains the NNLO mixed coulombic-logarithmic correction. Finally, 



/\(2) _ J2) PWs 



c^^\ ~ 1.75. 



(4.6c) 



For the exact value of c^^\, see the Appendix, Eqs. (A9). This finishes the calculation of the decay rate. 
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4.2 NLO calculation 

Using (4.1-6) to one loop only we find the NLO result, 



1 



3/3o I , , 
— ;- log ^ - 7E + 



r{T^e+e-) = 

xr(°)(r ^ e+e") 
where the LO expression F^^^ is 

r(o) = 2 



7e/3o ^ 93-10n/^^ 



36 



-2Cf 



(4.7) 



QbCtQED 



M{T) 



[mCFas{ij)f 



We note that (4.7) is slightly different from the result obtained using a variational method (instead 
of Rayleigh-Schrodinger perturbation theory) to evaluate 5R\o, which gives 



3/3o A «M \ , 
^(log--7Ej + 



7e/3o ^ 93-lOnj 



36 



2Cb 



^lr(°)(r^e+e-). 

TT 



4.3 NNLO calculation, and numerical results 

The NNLO result is obtained with the full (4.3) (4.6) expressions, with which numerical results are readily 
obtained. We select the values of the basic parameters A{nf = 4, three loops), {asG^) as before. For the 
renormalization point there are two "natural" choices: 

Hi = mb and = 2/a = 2.866 GeV . 

The first value gives a reasonable value for the width. These are the input values taken for the calculations 
reported in the Table 2 below. 



r{r e+e-) 


LO 


NLO 


NNLO 


IJ. = m 
fj, = 2/a 


0.41 
0.73 


1.22 
0.55 


1.13 kcV 
0.80 keV 



Table 2: Determination of r{T —f e+e ) to increasing accuracy. 



The experimental figure is 

r{T e+e-) = 1.32 ± 0.04 keV . 

Clearly, the large NLO and NNLO perturbative corrections, both of similar size, and of the leading NP 
correction, make the theoretical result unstable, as the Table shows. 

For the (perhaps measurable) toponium (T) width we get, for rut = 175 GeV, the corresponding 
results, summarized in Table 3. 





e+e-) 


LO 


NLO 


NNLO 


M = 


mt 


6.86 


10.53 


13.0 


M = 


2/a 


10.24 


10.91 


13.5 



Table 3: Determination of r{T e+e ) to increasing accuracy. 



Tlio situation has improved with respect to what we had for bottomium, but the dependence on the 
renormalization point and on the order of perturbation theory considered is still a bit large. We can however 
conclude on an estimate of some 11.0 — 14 keV for the width. 
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Appendix: Perturbation by terms p log" p 

To obtain the result of perturbations by terms p^^ log"p we first evaluate the perturbation by a term p'', 
V = integer, continue to noninteger v and differentiate with respect to v dX v = —1. For this we use the 
method devised by Leutwyler!''! and extended in the second paper of ref 5, Appendix, in general. The basic 
formula is the inversion formula 

here 

h{k) = Ho + RM = -Ir^^-'i" 

r a"*/^ 

with Hq the free hamiltonian and a = 2/mCpas. This inversion formula is valid for any k and integer p. 
Using this one may then check that the first order perturbation by p" , ^(y)-Rio, is 

, (Al) 



^+V(z^ + 2)-^(l) 



Sums like that in (Al) may be continued to arbitrary u by using the formula 

z^+l oo 

E/(^') = E + ^ + (A2) 

i=i j=i 

This replacement is valid provided the sum X^j^i fU) convergent; otherwise, we have to separate from / 
the leading, next to leading . . . pieces, to be summed explicitely, leaving a residue for which the sum up to 
infinity converges. This has already been done in getting (Al). As a check of this, we mention that the result 
agrees, for = —1, with the result of the (trivial) evaluation obtained directly by replacing the coulombic 
potential according to 

nCpots nCpots 1 
r r p 

and expanding: (5(_i)i?io(p) = — (ma^/4)(| — ip)i?io(/9). The results of first order perturbation by p~^ logp, 
p~^ log^ p are then found by differentiating (Al) with respect to v. 
To second order we require 



SliX.)RlO = PlO ^(0) _ ^(0) ^lo/PlO ^(0) _ ^^^^ PioP'Rio 



and 



^2{Xv)P'W — {Rw\p^\Rw) PlO — (0) ITTTr-'^iOTTfVi — ZZ7i:^^Op'^RlO- 
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Here Pio is the projector Pio = 1 — |i?io)(i?io|- At p = we then have, 

,„2^ 2 



<^l(Ai/)-^10(0) = - 2 

xi 3 



4 

r(A + 3) 



+ 1 



+ ^(1^ + 2) -^(1) 



r(4 + A) 



+2r(3 + A) 

,r(A + 5) 



1 



+ V(!^ + 2)-V(l) 



A + 1 



+ V(A + 2)-V(l) 



2 A + 2 

and the functions are defined as 



r{x + 4) 



A + l 



+ ^(A + 2)-V(l) 



^ jT(i + 2) 



jI^(A, ^) = ^ ^jT("j +tf [J' + ^ + 1 + ^ (^0' + ^ + 2) - ^(1))] 



J 



E ^ • nt-^ ^ [A + 1 + 2 (V>(j + A + 2) - ^(1))] 



jT(i + 2) 

r(A + ;/ + 5) 



r(A + 2)-r(A + 3), 



{X + 2)r{u + 3) 

for !/, A = integer. The second expression for ^ is obtained by using the identityl^^l 



E 

fc=0 



n + A: 



n + /X + 1 
n + 1 



and has the advantage that it can be continued directly with the use of (A2). 
We will require the following results: 

dvdX ^ ' A=i/=-i ^ 7^ 



<3 + l 



dvdX 



^{X,v) 



= 2-7E 



OC 

E {2i [^"(i + 1) + ^'(i + + 2)] + 2 + 2) - VO' + 2)) [^{3 + 1) - 

-2/(j + l)-l}//. 



With the same methods we find, 



5\x,^Rm = -\r{u + 3)r{X + Z) 



ma 



+ V(z^ + 2)- V(l) 



(A3a) 



(A3b) 



(A4) 



(A5) 



and we will need d^lQ.v)/ dv\y=-\. Write 



-15 + >f(0,z.)|i?io(0), 
(A6) 



^{Q,u) = -3 + 



{v + A){v + Z) 



+ v + l + 2[i,{v + 2)- V(l)] + 2 ^ + 2) - V(l)] + 4 ^ 



V(j + 2) - V(l) 



(A7a) 
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With (cf. ref. 11 for some of the harmonic- type sums here^ ) 



d 



-^[V;(j + 2)-^(l)] 



i=i 



we find 



d — V;(j+2)-V>(l) 

d^{0, v) 



^.=-1 3 



=C3 + 1 



(A7b) 



dv 



4 (Ca + 1) + TT^ - 1/2. 



(A7c) 



The function 

-Rio = -Rio + <^(i/)-RlO + ^\\v)Rw 

is not normaUzed to unity. The normaUzed wave function is i?io = ||-Rio||~^^^-Rio- Because i?io is orthogonal 
to 5i,R\o, ^\v^w the normaUzation factor is 

ll^ioir' = {1 + (5(.)i?lo|^(A)i?lo)}"'^' ^ 1 - W^.)Rw\SiX)Rw)- 
Therefore, we have to correct 5^^^^iiio for this and define 



"A+1 



r(A + 3)r(j^ + 3) fma 



X ! + i 



A + 1 



2 +V(A + 2)-V(l) + ^+V(i^ + 2)-V(l) 



■^(A + 2)-V(l) 



+ V(!^ + 2)-V(l) 



(A8) 



(z/+l)(z/ + 2) + (A + l)(A + 2) 



+ 5(z/ + A-2) 



i^+l A+1 



+6(V;(;/ + 2) + ^|J{X + 2) - 2^(1))] + ^ ^ ^ 



r(j + fc + 3) 

,tt^jAi+2)fcnfc+2) 

The only difiicult point is the continuation of the double sum 

i^+l A+1 



>R 



10- 



r(j + fc + 3) 

tfoi^(i + 2)fcr(fc + 2) 

which diverges if using directly (A2). We avoid this by resorting again to the identity (A4). With it, we get 



5«=V;(. + 3) + ^(A + 3)-2^(A + 2)-2-2[|^^^pi-^ + |^'l + l]+(:. + 3) (2 + J) 

fe— 2 fe— 2 



A+1 



- 1 y 



2r{v + fc + 2) _^ 2r(i/ + fc + 3) _^ 2r{v + fc + 4) 



fc 



fc + 1 



1 



fcr(fc + 2) ' 



fc- 1 

which can now be continued with (A2). 

^ The sum X^J^i ['AO + 2) — V'(l)] can be evaluated exactly using the identity (A4) differentiated with respect to n 

at n = 0. The derivative with respect to v of the sum Y^'^j'^l V'O + 2) — V^W ^ being only logarithmically divergent, 

3 

can be calculated with (A2). 
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For our calculation, we need to evaluate 



dvd\ 



v=\=-l 

The only nontrivial piece is the apparently divergent one arising for fc = 2. To get it, we note that 

92 1 2r{v + A + fc + 2) 



dvdX r{v + 3) (A + k){\ + k- l)r{X + k + 2 



|V'(2)-i<(2) = 7rV2 + 2C3-5. 



fe=2.7y=A=-l 



It is possible also to check some of the calculations here in two respects. First, we may evaluate the 
second order correstion to the energy shift by Hil ■ This is given by 



6rE,o = - 



2Y Capiat 



-{Rw,p^Pw- 



Using our formulas here we find 



-m 



47r2 diydX ^ E(O) - 
C%l3lal 3 + 37I - TT^ + 6C3 



PwP^Rw) 



47r2 12 

in agreement with the result of ref. 5, obtained with a completely diferent method (Green's functions). 
Secondly, we can particilarize the calculation here iov v = \ — —1. We get 



^f-i,-i)^io = (! - fp + \p') (^) i?io - IR. 



and the last term is the normalization correction. This is in agreement with the direct result 

-,2x 2 



'^^i,-i)^io = (i-|p+|p^)(^) • 



With these evaluations we can get the result of the second order perturbation by p ^ log p. Write 



g2 /^„2X 2 



.R, 



Rio [di 



One defines 
92 

No =^wT^Sx. 



dXdv 

00 { 

+2^ 2 
fe=i I 



\=v=-l 



+ 2C3 



V;(fc + 2) - V>(2) 
k{k + l){k + 2) 



1 



1 



1 



k k+1 k+2 



M + 2)-i^{2) 

fc3 



7^(fc + 3)-VX2) 
k{k + l) 

~10.29; 
92 



k + 2 



fc + l 



-2 



^'{k + 2) 
k{k+l){k + 2) 



,V'(fc + 2) (fc + 2)V''(A: + 3) 



A;2 



k{k + l) 



'dXdv 



^{X,u) 



6.82 



(the exact value of the last given in (A5)). Moreover, 



' \=u=-l ^ 



dXdu ^ ' ' ^ j2 ' 5 

Some of the sums in the numbers Ni are simplified in terms of the following harmonic sums'^^l: 



(A9a) 



j=l •' k=\ j=l •' fe=l 

00 00 ^ 

C3; E^E(^:p7)2='^V 



00 ^ 00 



120. 
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Then, 



rfi=-i|(l-7E) ^-7r2 + 67E- 
d2=- i(l - 7e) {-f + 67E + 4C3} ; 

dnorm. = " |(1 " 7e)' + f (1 " 7e) + M ^ + y 

Finally, the coefficient Cj^ \ is 



- 3 



C3 + I- 



7r^7E" 



-LA 



di- d2 + d 

n 



1.75. 



(A9b) 



(A9c) 
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